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Abstract. 



o 
o 

' We present a direct A-body simulation modeling the evolution of the old (7 Gyr) open cluster 

NGC 188. This is the first A?-body open cluster simulation whose initial binary population 
is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare 
the simulated color-magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue 
stragglers produced in the simulation. We compare the solar-type main sequence binary period 
• and eccentricity distributions of the simulation to detailed observations of similar binaries in 
£\j | NGC 188. These results demonstrate the importance of detailed observations in guiding A-body 
open cluster simulations. Finally, we discuss the implications of a few discrepancies between the 
NGC 188 model and observations and suggest a few methods for bringing Af-body open cluster 
Ph ' simulations into better agreement with observations. 
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■ 1. Introduction 

Recently, sophisticated A^-body simulations have incorporated stellar dynamics and 
stellar evolution self consistently, gaining the ability to accurately model open cluster 
-H ■ sized populat ions (N ~ 10 ) of single and binary stars through many bill ions of years . 
(e.g. NB0DY 6, Aarsethll2003l with stellar and binary evolution included by iHurlev et al 



20001 l2002h . Concurrently, large observation al surveys of open clusters, like the WIYN 
Open Cluster Study (WOCS: lMathieull2000l ), are coming to maturity, providing compre- 
hensive databases of open cluster characteristics, including detailed information on their 
binary populations. These observations, and specifically those of the binaries, provide 
important tests and guidance for A^-body simulations. In this paper, we describe our 
^\ ' first step in utilizing this wealth of data to help guide an A^-body simulation with the 
■ aim of creating a realistic model of NGC 188. 

2. TV-body Model of NGC 188 



The groundbreaking IHurlev et al. (2005) model of M67 is a clear example of the po- 



tential for Af-body simulations to model real open clusters, and thereby address long- 
standing questions about the evolution of a binary population, origins of blue stragglers 
(BSs), and specific observed characteristics of a given open cluster. Hurley et al. had very 
limited knowledge of the real binary population in M67, and thus aimed to match other 
important observables like the number of BSs, the spatial profile, etc. Today through 
the WOCS su rvey, the work at the Harvard-Smithsonian Center for Astrophysics (e.g. 



Latham 2006) and others, we now know a great deal about the M67 binaries. It is 
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Figure 1. Period (left) and eccentricity (right) binary frequency distributions comparing the 
primordial binary population in the simulation (filled histograms) to the observed MS binary 
population of M35 (150 Myr; solid black lines). The M35 observed distributions have been 
corrected for incompleteness to show all binaries with P < 10 4 days (excluding those with 
periods below the circularization period in the eccentricity distribution plot). The primordial 
binary distributions used in the simulation are plotted with dark gray for observationally de- 
tectable binaries (P < 10 4 days), and light gray for longer-period binaries. We also show the 
Duqucnnoy & Mayor (1991) Galactic field binary frequency distributions in the dashed lines for 
comparison. 



now clear that, though iHurlev et al. (2005) were successful in matching many observed 
characteristics of M67, the binary frequency at 4 Gyr was too high and the simulation 
contained too many sh ort-period binaries throughout the evolution of the cluster (see 
Geller fc Mathieul l2009) . The binary population has a significant effect on the formation 
mechanisms and rates for BSs as well as the dynamical evolution of the entire clus- 
ter. Therefore, correctly modeling the binary population is one of the most important 
ingredients in creating a realistic open cluster model. 

The old (7 Gyr) open cluster NGC 188, like M67 (4 Gyr), has a r ich and well-studied 
binary population (e.g., Geller et al. 2009a : Geller fc Mathieu 2009), a large population 
of BSs dMathieu fc Gellerl l2009h. " and a wealth of observations from the WOCS group 
and others (see lGeller et al.ll2008l ). Indeed, NGC 188 is quickly becoming one of the most 
well studied and important open clusters for understanding the origins of BSs, the char- 
acteristics of an evolved binary population and the dynamical state of an evolved open 
cluster. These data provide ideal and much-needed guidance for an iV-body simulation 
of NGC 188. 



2.1. Simulation Method 

All simulations were run using the NB0DY6 code (|Aarsethl l2000). with slight modifications 
added by the authors to define the initial binary po pulation and output format. All stellar 
evolution is incorporated as in Hurley et al. (l2000f) and b inary evolution, including tides, 



mass transfer 
manner as in 



etc., is trea t ed as in IHurlev et al. (2002). BSs are modeled in the same 



Hurley et al.1 (|2005l ). The simulations were executed on the supercomputer 
at the Swinburne Centre for Astrophysics and Supercomputing. 

The NGC 188 simul ation was initi alized with a mass of 23,870 M© (N = 24,000 stars) 
chosen from an IMF of lKroupal (]200 lh and distribute d according to a Plum mer profile. We 
assume a solar metallicity for the stars in the model ( Saraiedini et al.lll999h . and place the 
cluster in or bit around a point-mass Galactic potential according to the observed orbit 
of NGC 188 (jCarraro fc Chiosilll994f) . In the following, we discuss the binary population 
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of the model in some detail, focusing on the period and eccentricity distributions and 
leaving more thorough discussions of methods and results for future papers. 
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Figure 2. Color-magnitude diagrams of NGC 188 (left) and the NGC 188 simulation at 7 Gyr 
(right). The single stars are plotted in the small black points, and the detected binaries are 
plotted with the l arger open circles. The NGC 188 BSs are identified as being to the left of the 
dashed line (as in lGeller et al.|[2008f ). and within the simulation BSs are identified as being MS 
stars with masses greater than the turnoff mass (and reside in a similar region to observed BSs 
in NGC 188). The simulation succeeds in reproducing the overall form of the observed single 
and binary star sequence, but has not created the large number of observed BSs. 



2.2. Primordial Binary Population and Observations of M35 

Choosing the initial parameters defining a binary population has long been derived from 
theoretical analyses and predictions, as we have lacked sufficiently detailed observations 
of young open cluster binary populations to help guide these choices. The WOCS survey 



of M35 (|Geller et alJ l2009bh now allows for a major step forward; we can use these 



binaries as a proxy for the primordial binary population in this simulation. There is no 
a priori reason to presume that the M35 binaries will evolve into the NGC 188 binaries. 
Still, the detailed observations of the M35 binaries allow us to create the first TV-body 
model whose initial binary population is directly defined by observations and to study 
the evolution of these binaries for the entire age of NGC 188 (7 Gyr). 

M35 is observed to have a main-sequen ce (MS) binary frequency of 22 ± 3 % out 



to periods of 10 4 days (jGeller et al. 2009b); we have chosen to initialize the simulation 



with a binary frequency of 25% out to the same period range. In Figure [T] we show the 
M35 incompleteness corrected binary period and eccentricity distributions (black line) 
as compared to the initial conditions used in the NGC 188 model (filled histograms). As 
we only detect binaries wit h periods < 10 4 days, we utilize the Galactic field binaries 
(|Duquennov fe Mavorlll99ll ) to extend the period distribution to longer periods. Thus 
we input a primordial binary population with a log-normal period distribution and a 
Gaussian eccentricity distribution. 

2.3. Comparing to Observations of NGC 188 

The WOCS observations of NGC 188 (|Geller et al.ll2008ll2009aHGeller fc Mathieull2009h 
allow for a detailed comparison with the NGC 188 model. In the following analysis, we 
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Figure 3. Period (left) and eccentricity (right) MS binary frequency distributions at 7 Gyr. The 
NGC 188 observed distributions are plotted in the solid black lines and have been corrected for 
incompleteness in the same manner as in Figure [T] showing all binaries with P < 10 4 days and 
excluding those with periods below the circularization period in the eccentricity distribution 
plot. The simulated binary distributions are plotted in the filled histograms with the same 
gray-scale coding as in Figure [1] Again, we also show the iDuquennov fc Mayor! l|1991f) Galactic 
field binary frequency distributions in the dashed lines for comparison. Note the close agreement 
between the simulated and observed binaries in both distributions and binary frequency. 



limit the output of our simulation to only include stars within the Geller et al. observable 
domain (in both mass and binary period) . We first focus on the color-magnitude diagram 
(CMD) at 7 Gyr shown in Figure [U Note that the stellar evolution incorporated into 
the simulation does an excellent job of reproducing the "normal" single and binary stars 
sequence as compared to the observed CMD. However, the simulation is notably lacking 
in BSs; there are 21 observed BSs in NGC 188 and only 5 at 7 Gyr in the simulation. 
This low number of BSs stays relatively constant throughout the entire simulation, and is 
likely due to the initial binary population. We will discuss the implications of this result 
in Section [31 

We find the solar-type MS binary population in the simulation at 7 Gyr to correspond 
remarkably well to observations of similar binaries in NGC 188. In Figure [3]we compare 
the observed incompleteness-corrected NGC 188 binaries (black lines) to those of the 
simulation at 7 Gyr (filled histograms) in both period and eccentricity. The form of 
both the simulated distributions, as well as the binary frequenc y agree well with the 
obser vations of NGC 188 as well as the Galactic field (dashed lines; IDuquennov fc Mayor 



19911) . This close agreement is a testament to the power of using detailed observations 



to define the initial binary population in a A^-body simulation. 

Still, there are a few discrepancies between the simulated binary population and the 
observations. The circularization perio d in the model is si g nifica ntly shorter than the 
observed value for NGC 188 (Figure l4j). iMeibom fc Mathieul (|2005f ) find a circularization 
period of 14.5+2 2 days, while a similar analysis of the model shows a period of 2.5 days. 
Interestingly, we also see too many circular MS binaries in the simulation at 7 Gyr (first 
bin in Figured]), most having periods greater than the circularization period of both the 
simulation and observations. We discuss both of these results in the following section. 



3. Discussion and Future Directions 

By utilizing the observed M35 binaries as a proxy for our primordial binary popula- 
tion, we have managed to nearly reproduce the NGC 188 binaries within our A^-body 
simulation at 7 Gyr. However, as noted in Section [2~3l the current NGC 188 model is de- 
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Figure 4. Eccentricity plotted against the logarithm of the period (e - log P) at 7 Gyr for 
the observed (top) and simulated (bottom) solar-type M S NGC 188 binaries . The solid gray 
lines represent the best-fit circularization functions from [jylcibom & Mathicu, ()2005l ) for each 
sample, respectively. The tidal processes in the simulation result in a circularization period that 
is significantly lower than in NGC 188. Also, the simulation produces a population of circular 
binaries with P » P C i rc that is not observed in NGC 188 (or the Galactic field). Most of these 
circular binaries in the model have previously gone through a stage of mass transfer. 



ficient in BSs. Interestingly, the iHurlev et al] (|2005h model of M67 was able to produce 
the large number of observed BSs (though not their binary frequency), but was unable 
to reproduce the observed binary population. The quantity of BSs produced in the M67 
model was largely dependant on the initial binary population. Hurley et al. found that 
a large fraction of short-period binaries was necessary to create the observed numbers of 
BSs in M67 (and NGC 188). These short-period binaries not only have a higher likeli- 
hood of creating BSs in isolation (e.g., through mass transfer resulting from Roche lobe 
overflow), but dynamical interactions involvi ng short-period binar ies have a higher prob- 
ability of resulting in collisions and mergers (jFregeau et al.ll20041 ). Observations of M35 
binaries are very different from the flat period distribution and large binary frequency 
used in the M67 model (e.g., Figure QJ. The important next step will be to match both 
the binary population and the BSs population within one cohesive simulation. We are 
currently exploring a number of possible methods for addressing this issue, the most 
promising of which appears to be the inclusion of primordial triples (e.g., Kozai induced 
mergers; Perets fc Fabrvckv 20091 . larger cross section for dynamical interactions, etc.). 

The tidal circularization period is one of our best observational tools to study the 
effects of tides on a binary population. In Figure [Hand Section |2~51 we show that the cir- 
cularization period at 7 Gyr in the simulation is significantly lower than what is observed 
in NGC 188. This suggests that the tidal energy dissipation rate may be underestimated 
in the iV-body model (at least for solar-type MS binari es). We may also (or alterna- 
tiv ely) require a fo rm of pre-MS tidal circularization (e.g. Kroupal[l995l ). As pointed out 
bv iMathieul (|2008t ), tidal processes have a significant impact on dynamical interactions 
between single and binary stars, and the energy dissipation rates can be the deciding 
factor between a close fly-by, tidal capture, or even a merger or collision that could result 
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in a BS. WOCS data (e.g., Meibom fc Mathieull2005 ) provide the ideal foundation from 



which we will improve the tidal physics within the TV-body model. 

Finally, the simulation has an excess of circular binaries with periods beyond the cir- 
cularization period. The majority of these binaries have undergone some form of mass 
transfer, which has quickly circularized the orbit within the simulation. However, similar 
populations of circular binaries at P » P C j rc are not observed in open clu sters or in 
the Galactic field (e.g.. lMeibom fe Mathieul 2005: Du quennov fe Mavorlll99l[) . This may 



suggest that mass transfer, e specially in initially eccentric binaries, does not necessarily 
lead to circular binaries (e.g., Sepinskv et al. 2007 : Bonacic Marinovic. Glebbeek fe Polsl 



2008f ). or that these cases of mass transfer happen much more frequently within the sim- 



ulation than in reality. Correctly modeling the products of mass transfer and common- 
envelope evolution is critical for accurately reproducing BSs as well as the binary popu- 
lation within an TV-body simulation. 

We will present further details about this NGC 188 model in future papers. Addition- 
ally, we will explore the various methods mentioned above for bringing the NGC 188 
model, and iV-body simulations in general, in better agreement with observations. This 
synergy between observations and simulations will continue to refine the TV-body method 
and reveal further insights into the origins of BSs, the evolution of a binary population, 
and the dynamical evolution of star clusters. 

This work was funded in part by a National Science Foundation (NSF) East Asia and Pacific 
Summer Institutes (EAPSI) fellowship, the Wisconsin Space Grant Consortium, and NSF 
grant AST-0406615. 
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